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1 Introduction 

Fractional Brownian Motion (F.B.M.) was introduced in 1940 by Kolmogorov as a way to generate Gaussian 
"spirals" in a Hilbert space. But the seminal paper of Mandelbrot and Van Ness (1968) emphasizes the rel- 
evance of F.B.M. to model natural phenomena: hydrology, finance... Formally, a fractional Brownian motion 
Bh = (Bu(t), t 6 could be defined as a real centered Gaussian process with stationary increments such that 
B H (Q) =0and E \B H (s) - B H {t)\ 2 = a 2 \t-s\ 2H , for all pair (s,t) G M + x M + where H G]0, 1[ and a > 0. This 
process is characterized by two parameters : the Hurst index H and the scale parameter a. We lay the emphasis 
on the fact that the same parameter H is linked to different properties of the F.B.M. as the smoothness of the 
sample paths, the long range dependence of its increments and the self-similarity. 

During the decades 1970's and 1980's, the statistical study of F.B.M. was developed, to look at for instance 
the historical notes in Samorodnitsky & Taqqu (1994), [301 chap. 14] and the references therein. Modelling by a 
F.B.M. became more and more widespread during the last decade (traffic Internet, turbulence, image processing...). 
Nevertheless, in many applications the real data does not fit exactly F.B.M. Thus, the F.B.M. must be regarded 
only as an ideal mathematical model. Therefore, various generalizations of F.B.M. have been proposed these last 
years to fill the gap between the mathematical modelling and real data. In one hand, Gaussian processes where 
the Hurst parameter H has been replaced by a function depending on the time were studied, see for instance 
Peltier and Levy Vehel (1996), Benassi, Jaffard and Roux (1997), Ayache and Levy Vehel (1999). However, this 
dependence of time implies the loss of the stationarity of the increments. In other hand, non Gaussian processes, 
mainly a stable (0 < a < 2) infinite variance processes, were considered, see for example the study of telecom 
processes in Pipiras and Taqqu (2002). 

Here, we are concerned with Gaussian processes having stationary increments and a Hurst index changing with 
the frequencies. To our knowledge, these kinds of processes were introduced implicitly in biomechanics by Collins 
and de Luca (1993), in finance by Rogers (1997) and Cheridito (2003) and explicitly by Benassi and Deguy (1999) 
for image analysis or image synthesis. In any case, the probabilistic properties of these processes have not been 
thoroughly established and no rigorous statistical studies have been done. Both Collins and de Luca (1993) and 
Benassi and Deguy (1999) propose a model with two different Hurst indices corresponding respectively to the high 
and the low frequencies separated by one change point at the frequency u> c . They use the log variogram to estimate 
these two Hurst indices. Indeed, in this case, the log variogram considered as a function of the logarithm of the scale 
presents two asymptotic directions with slopes being twice the Hurst index at low (respectively high) frequencies. 
The change point u> c is then estimated as the abscise of the intersection of the these two straight lines. Numerically, 
this method is not robust. Moreover it could not be adapted in the case of more than one change point. Let us 
stress that it is not a question of a theoretical refinement, but one that corresponds precisely to the true situations. 
Indeed, in applications, we consider only finite frequency bands, therefore we should use a statistical method based 
on the information included in finite frequency bands. Wavelet analysis seems the tool had hoc, when the Fourier 
transform of the associated wavelet is compactly supported. 

For these reasons, we put forward in Bardet and Bertrand (2003) a model of generalized F.B.M. including the 
cases with more than one frequency change point. We called it (Mr-) multiscale fractional Brownian motion where 
K denote the number of frequency change points. More precisely, a (Mr) multiscale fractional Brownian motion 
is a Gaussian process with stationary increments where the Hurst parameter H is replaced by a piecewise constant 
function of the frequency £ <— > H(£) in the harmonizablc representation, see Formula (J3]) below. The main proba- 
bilistic properties of this model were studied in Bardet and Bertrand (2003). In this work, we treat the statistical 
study of the multiscale F.B.M. and we focus on its application to biomechanics. 
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The remainder of the paper is organized as follows: in Section 2, we describe the biomechanical data and the 
corresponding statistical problem. In section 3, we recall the initial definition of the partial Brownian motion and 
its principal probabilistic properties. Then, we show that the variogram method is not suitable for the estimation 
of the various parameters of a (Mk)-F.B.M. We then develop a statistical estimation framework, based on wavelet 
analysis. We investigate the discretization of the wavelet coefficient and we state a functional Central Limit The- 
orem for the empirical wavelet coefficients. In Section 4, we first estimate the different frequency change points 
and Hurst parameters. Then, we propose a goodness of fit test and derive an estimator of the number of frequency 
changes. The numerical algorithm is detailed at the end of this section. Finally, in Section 5, the biomechanical 
data are studied with the tools developed in Section 4. The proof of the results of Sections 3 and 4 are given in 
appendix. 

2 The Biomechanical Problem 

One of the motivations of this work is to model biomechanical data corresponding to the regulation of the upright 
position of the human being. By using a force platform, the position of the center of pressure (CO. P.) during 
quiet postural stance is determined. This position is usually measured at a frequency of 100 Hz for the one minute 
period, which yields a data set of 6000 observations. The experimental conditions are formed to the standards of 
the Association Francaise de Posturologie (AFP), for instance the feet position (angle and clearance), the open or 
closed eyes. 



Figure 1 : An example Q| of the trajectory of the CO. P. during 60s at 100Hz (in mm) 

The X axis of the platform corresponds to the fore-aft direction and the Y axis corresponds to the medio-lateral 
direction. During the 1970's, these data were analyzed as a set of points, i.e. without taking into account their 
temporal order. During the following decade some studies considered them as a process, and Collins and de Luca 
(1993) introduced the use of F.B.M. to model these data. In fact, they used a generalization of F.B.M. More 
precisely, let the position Xi of the CO. P. be observed at times i, = iA for i = 1, . . . , N (A = 0.01 s). The study 
of Collins and de Luca is based on the empirical variogram 
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Position following X axis 




N-S 



(1) 



1 these experimental data were realized by A. Mouzat and are used in |13| . 
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where 8 e M*. For a F.B.M., we have EVn(S) = a 2 A 2H x S 2H and after plotting the log-log graph of the 
variogram as a function of the time lag , i.e. (log 8, log Vn (#)), a linear regression provides the slope 2H. Typically, 
one gets the following type of figure (see Figure 2). It is considered by Collins and de Luca to be a "F.B.M." with 
two regimes : with slope 2Hq (short term) and with slope 2Hi (long term) separated by a critical time lag S c and 
these parameters are estimated graphically : 




_4 1 1 1 1 1 1 1 1 

1 2 3 4 5 6 7 8 

Logarithms of the scales 

Figure 2 : An example of the log-log graph of the variogram for the previous trajectories X (-.) and Y (-). 

They found Hq > 0.5, H\ < 0.5 and a critical time lag 5 C ~ 1 s. These results were interpreted as corre- 
sponding to two different kinds of regulation of the human stance : in the long term Hi < 0.5 and the process is 
anti-persistent, in the short term Hq > 0.5 and the process is persistent. This method was employed several times 
in biomechanics under the various experimental conditions (opened eyes versus closed eyes, different feet angles,...). 
But, a lack of mathematical models and of statistical studies has made impossible to obtain confidence intervals 
on the two slopes 2H , 2H\ and the critical time lag S c . 

3 The multiscale fractional Brownian motion and its statistical study 
based on wavelet analysis 

3.1 Description of the model 

A fractional Brownian motion Bh = {Bn(t), t e M} of parameters (H, a) is a real centered Gaussian process with 
stationary increments and E \B H (s) - B H (t)\ 2 = a 2 \t - s\ 2H , for all (s,t) € M 2 where H e]0, 1[ and a > 0. The 
fractional Brownian motion (F.B.M.) has been proposed by Kolmogorov (1940) who defined it by the harmonizablc 
representation : 

B « (*) = J R \\H+i/2 fOT a11 1 e K ' ( 2 ) 

where W(dx) is a Brownian measure and W{d^) its Fourier transform (namely for any function / G L 2 (M) 
one has almost surely, J R f(x)W(dx) — j M /(£) W(dt;), with the convention that /(£) = Jj R e" l ^ x f(x) dx when 
/ e L 1 (M) P| L 2 (M)) . We refer to Samorodnitsky and Taqqu (1994) for the question of the equivalence of the 
different representations of the F.B.M. From the harmonizable representation, a natural generalization is the 
multiscale fractional Brownian motion with a Hurst index depending on the frequency. More precisely, we define : 
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Definition 3.1 For K S IN, a (MK)-multiscale fractional Brownian motion X = {X(i),t £ M} (simplify by 
(Mk)-F.B.M.) is a process such as 

X(t) = 2j2 *j \ fl H ]+ i/2 W ( d for all tefl (3) 

with luq = < uj± < ■ ■ ■ < ll>k < Wjf+i = oo by convention, c*i > and Hi £E]0, 1[ for i G {0, 1, • • • , K}. 

The (Mr-)-F.B.M. was notably introduced in order to relax the self-similarity property of F.B.M. Indeed, the 
self-similarity is a form of invariance with respect to changes of time scale [37] and it links the behavior to the high 
frequencies with the behavior to the low frequencies. In Bardet and Bertrand (2003), the main properties of these 
processes are provided : X is a Gaussian centered process with stationary increments, its trajectories are a.s. of 
Holder regularity a, for every < a < Hk and its increments form a long-memory process (except if the different 
parameters satisfy a particular relationship, i.e., if its spectral density is a continuous function with < H \ < 1/2 
for i = 0,1,- • • ,K). 

3.2 The question of the choice of the estimator 

In the remainder of this paper, we suggest a statistical study of such a model based on wavelet analysis. In this 
subsection, we explain the reason of this choice. 

To begin with, we will describe the statistical framework precisely. Let X = {X(t),t € be a (M^)-F.B.M. 

defined by Q. We observe one path of the process X on the interval [0, Tjv] at the discrete times U — i ■ Ajy for 
i=l,...,N with Tjv = N ■ A N . Therefore, 

(X(A N ),X(2A N ), . . .,X(NA N )) is known, 

and we consider the asymptotic N — » oo, An — * and T/v — > oo. We want to estimate the parameters of the 
(M K )-F.B.M. that arc {H ,H l7 . . .,H K ), (a 0) ai, . . .,<tk) and . . .,uj k ). 

Even if the model is defined as a parametric one, we prefer to use a semi-parametric statistics based on the 
wavelet analysis. This choice is justified by the following reasons. First, the spectral density of X is not continuous 
in the general case. Thus, one cannot use the classical results on the consistency of the maximum likelihood or 
Whittle maximum likelihood estimators for long memory processes (see Fox and Taqqu, 1986, Dahlhaus, 1989 or 
Giraitis and Surgailis, 1990). Moreover, this is not a classical time series parametric estimation : indeed, we con- 
sider (X(A N ),X(2A N ), . . .,X(NA N )) instead of (X(1),X(2), . . .,X(N)) and therefore this is also an estimation 
problem of the parameters of a continuous stochastic process. Secondly, the following semi-parametric statistics are 
more robust than a parametric one if the model is misspecified. Consider the example where the function H (£) is a 
not exactly a piece- wise constant function, but instead a constant function on several intervals and some unknown 
function on the other intervals. In this case, a parametric estimator could not work while the semi-parametric 
method based on the wavelet analysis will remain efficient. 

Another semi-parametric method was developed from the seminal paper of Istas and Lang (1997). This method 
of estimation is derived from the variogram and provides good results in the case of F.B.M. (see Bardet, 2000) 
or of multifractional F.B.M. (see Benassi et al., 1998). However, one faces difficulties in identifying the model 
(Mif)-F.B.M. with this kind of method. Indeed, one can easily satisfy that for S > : 

V(S) = E(X(t + S)-X(t)f = 4j2 S2Hi °* / 2hZ dv - ^ 

3=0 JSuJ J V 
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The principle of the variogram's method ensues from the writing of log ( V(*)J as an affine function of log 6. For a 

(1 — cosv) 



(M K )-F.B.M., with C(Hi) = 



■ dv for i = 0, 1, . . . , K , two cases could provide such a relation : 



1. for 5 -» oo, log (V(*)) = 2H • log* + log (4 • a 2 ■ C(H )) + 0(<T 2H °); 



2. for 5^0, log V(*) = 2H K • log J + log (4 • a£ • C(#* )) + 0(<S 2 - 2H * ) 



(the proof of such expansions is in the proof of Lemma lA.l|) . In those cases, if one can show that there is a conver- 
gent estimator Vn(S) of V(<5), then a log- log regression of log \ Vn(S) \ onto log* could provide an estimation of the 
different parameters. Nevertheless, such a method would have a lot of drawbacks. On one hand, the estimation 
of "intermediate" parameters (Hj)i<j<K-i and (<7 2 )i<j<.R--i requires very specific asymptotic properties between 
all the frequency changes (Wj)i<j<ic— i. This implies a lack of generality of the methods based on the variogram. 
Moreover, concretely, the frequency changes are fixed and one obtains rough approximation instead of asymptotic 
properties. For instance, numerical simulations show that in some cases the log-log plot of the variogram does not 
exhibit any intermediate linear part. On the other hand, when the model is misspecificd the variogram model could 
lead to inadequate results. For example the following picture gives the case of a (M2)-F.B.M. where the variogram 
method would detect only one frequency change and could not precisely estimate its value. Finally, the variogram's 
method could perhaps be applied in the two first previous situations 1. and 2., i.e. for the estimation of (i?o,cr 2 ) 
or (HktOk) w hh * will have to be a function of N (number of data). But this choice of function will depend on 
the unknown parameters Hq or Hk for obtaining central limit theorems for log (Vn(8)^--- (see the same kind of 
problem in Abry et at, 2002). 




o 

lOQ(d) 



Figure 3: An example of a theoretical variogram for a (Af2)-f-B.m, with Hq = 0.9, Hi = 0.2, H2 = 0.5, and 
(To = <7i = (T2 = 5 and u)\ — 0.05, uj% = 0.5 (in solid, the theoretical variogram, in dot-dashed, its theoretical 
asymptotes for * — * and * — * 00). 



We deduce from the definition of the model and the previous discussion that a wavelet analysis could be an 
interesting semi-parametric method for estimating the parameters of a (M^)-F.B.M. Indeed, such a method is 
based on the change of scales (or frequencies). Therefore, as it is developed below, a wavelet analysis is able to 
detect the different spectral domain of self-similarity and then estimate the different parameters of the model. 
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3.3 A statistical study based on wavelet analysis 

This method has been introduced by Flandrin (1992) and was developed by Abry et al. (2002) and Bardet et al. 
(2000). We also use in the following similar results on wavelet analysis for (Mjf)-F.B.M. obtained in Bardet and 
Bertrand (2003). Let ip be a wavelet satisfying the following assumption : 

Assumption (Al): ip : M>—>Misa C°° function satisfying : 

• forall meM, / \t m ip{t)\ dt < oo; 

Jr 

• its Fourier transform ip(Q is an even function compactly supported on [— fj, —a] U [a, 8] with < a < 8. 

We stress these conditions are sufficiently mild and are satisfied in particular by the Lemarie-Meyer "mother" 

wavelet. The admissibility property, i.e. / ip(t)dt = 0, is a consequence of the second one and more generally, for 

Jr 

all me IN, 

t m 4>{t)dt = 0. (5) 



IR 

Note that it is not necessary to choose ip to be a "mother" wavelet associated to a multiresolution analysis of 
Ei 2 (M) . The whole theory can be developed without resorting to this assumption. The choice of ip is then very 
large. 

Let (a, b) € M* + x M and denote A = (a, b). Then define the family of functions ip\ by ip\(t) = —7=ip ( b 

Parameters a and b are so-called the scale and the shift of the wavelet transform. Let us underline that we consider 
a continuous wavelet transform. Let dx{a,b) be the wavelet coefficient of the process X for the scale a and the 
shift b, with 

d x (a,b) = 4= / i>{-~b)X(t)dt =< i/j x ,X > L 2 (R) . 

If ip satisfies Assumption (Al) and X is a (Mk)-F.B.M., the family of wavelet coefficients verifies the following 
properties (see Bardet and Bertrand, 2003) : 

1. for a > 0, (dx(a, fe))bgK is a stationary centered Gaussian process such as : 

E (d 2 x (a,.)) =li(a) = a [ $(au)\ 2 ■ p~ 2 (u) du. (6) 

JR 

a 8 

2. for alH = 0, 1, ■ • • , K, if the scale a is such as [— , — ] C o;,,o;,+iJ, then 

a a 

^. 2 

r *P( U ) 

E (d 2 x (a, .)) = a 2H ' +1 ■ u\ ■ K H M, with K H (^) = / T1MTT du. (7) 

Jr. \ u \ 

Property {7} means that the logarithm of the variance of the wavelet coefficient is an affine function of the logarithm 
of the scale with slope 2Hi + 1 and intercept log erf + log-fO/i (V0- This property is the key tool for estimating the 
parameters of X. Indeed, if we consider a convergent estimator of log (E (d x (a, .))), it provides a linear model 
in log a and logo - ?. Before specifying such an estimator, let us stress that one only observes a discretized path 
(X(0),X(An), . . . , X(NAjy)) instead of a continuous-time path. 

As a consequence, for a > and N 6 IN*, a natural estimator is the logarithm of the empirical variance of the 
wavelet coefficient, that is log In (a) where : 

Mo) = [p \ ni d 2 x (a,kA N ), (8) 

1 N ^ keD N (a) 

with : 
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• re]0,l/3[; 

• ton = [r(N/d)] and Mn = [(1 — r)(N/a)) where [x] is the integer part of x <E M: 

• D N (a) = {m N ,m N + 1, . . . ,M N } and \D N (a)\ is the cardinal of the set D N (a). 

For < a min < a max , a functional central limit theorem for (log/Ar(a)) a7llirl < a < amax can be established (see a 
similar proof in Bardet and Bertrand, 2003) : 

Proposition 3.1 Let X be a (Mk)-F.B.M., < a m in < a max and ip satisfy Assumption (Al). Then : 
VNA N (log I N (a) - log Ji(a)) a . 

„<a<a max > {^{ a ))a m in<a<a max (9) 

- - AT-+00 



wii/i (-Z(a)) a centered Gaussian process such as for (01,02) € [a m i„,a 



(l-2r)J 1 (a 1 )J 1 (a 2 )i K W K |p(£)| 2 / V ' 



Then, if we specify the locations of the change points in terms of scales, i.e. frequencies, we obtain the following: 

Corollary 3.1 Let i G {0, 1, • • • , K} and assume that - < -!±i. Then, 

a 0Ji 



^/NAn~ (log I N (l/f)+(2H l+ l) log /-log of -log K Hi ty)) 



w i/"</<^i+l//3 



V 



{Z(l/f)) Ui/a <f< Ui+1 /P (11) 
JV^oo 

with the centered Gaussian process (Z(.)) such as for (/1, / 2 ) G [— , —^—] 2 ! 

For An small enough, this result shows that all parameters Hi and of could be estimated by using a linear 
regression of log/jv(l//j) versus log/,, when the frequencies u>i are known. Moreover, this central limit theorem 
shows that a graph of (log/, log/jv(l//)) for / > exhibits different areas of asymptotic linearity : it suggests 
the procedure of the following section to estimate and test the frequency changes (see for instance figures 4 or 6). 

3.4 The discretization problem 

In the applications, we only observe a finite time series (X(0), X(An), • • • , X((N— 1) x Ajv)) and we must derived 
the empirical wavelet coefficients from this time series. Since the process X has almost a continuous path but with 
a regularity ax < 1 almost surely, we should use the Riemann sum. Thus, for (a, b) G M* + x It we define the 
empirical wavelet coefficient by 

e x (a,6) = ^^^-6)xI(pA N ) (13) 

V n a 

and the discretized estimator by 

Jn{o) = Yl e 2 x (a,kA N ). (14) 

1 N{ n keD N (a) 

We also define for every k G Djy(a) the error 

£jv(a,fc) = e x (a,fcAAr) - d x (a,kA N ). (15) 

Now, it is possible to provide the functional central limit theorem for (log ^w(a))a mi „<a<a max computed from 
(X(0),X(A N ),..- ,X(NA N )): 
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Theorem 3.1 Under assumptions of Provosition \3.1\ and with An such as NAn oo and N(An) 2 — > when 
N — * oo. Then, with the same process Z than in (0J) 7 

> /]^(bgJ^(o)-logI 1 (o)). 4n<0< _ ^ (Z{a)) aml n<a<a max - (16) 



j4s a particular case, for i £ {0, 1, • • • , K\ and if — < , t/ien 



^/NA^ (log J JV (l//)+(2ff i +l) log/-loga?-log K H M)„ i/a < f < Ui+1 
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^(l//))^/^/^ (17) 

iV— >oo 



The convergence rate of the central limit theorem (|16|l is y/NApj- Thus, the discretization problem implies that 
the maximum convergence rate is o(N 1 ^ 4 ) from the previous conditions on An- 

4 Identification of the parameters 

First, let us describe the method on a heuristic level. From Proposition 13. 1[ Formula (jTTJ) , we have 

\ogJ N (l/f) = -(2H l + l) x log(/)+log(<7?) +log(K H M)+4 N) > ( 18 ) 
for the frequencies / which satisfy the condition 

log ( Ui ) - log(a) < log (/) < log (uj 1+1 ) - log(P). (19) 
Moreover we have (NAn) 1 ^ 2 (e*. J — -> fj))i<j<m- Formula (fig)) and condition p9[) mean that for 

\ ' i J l<j<m AT— toe ~ ~~ 

l°g(/) S [log ( w i) — l°g'( a )j log ( w i+i) — l°g(/3)]j we have a linear regression of logJjv(l//) onto log(/) with slope 

— {2Hi + 1) and intercept 

log of + logKHi{tp) and for log(/) £ [log (wj+i) — log(a), log (wi+2) — log(/3)] a linear regression with slope 

— (2Hi + i + 1) and intercept log of +1 + \ogKji i+1 (ip)- This is a problem of detection of abrupt change on the 
parameters of a linear regression, but with a transition zone for log(/) S ]log {u>i+i) — \og((3), log {wi+i) — log(a)[. 

Remark 4.1 Condition H9\) implies that Wj+i > — x Wj. Therefore we could only detect the frequency changes 

a 

sufficiently spaced. For instance, if we choose the Lemarie-Meyer wavelet, we get (3/a = 4 which leads to the 
condition uJi+\ > 4 x u>i. 

In this section, we describe the estimation of the parameters and a goodness of fit test. Both of them are based on 
the following assumption : 

Assumption (Bk) • The process X is a (M^-niultiscale fractional Brownian motion. This process is 
characterized by the parameters fl*, H* and cr* where ft* = (&;*,••• ,u* K ) with H* = (Hq , iJjf, . . . , H^) and 
cr* = (cr*,, <t*, . . . , c* K ). Moreover the following conditions are fulfilled 

• uf, -1 > — x for % = 1, • • • , K — 1; 

+ a 

• o^-t) { - h *y + k - <) 2 } > ° &nd 

• there exists a compact set JC c]0, l[x]0, 00 [ such as (if*, cr*) G /C for all i = 0, 1, • • ■ , K . 
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4.1 Estimation of the parameters 

Let X be a (Mk )-F.B.M. satisfying the assumption (Bk) with K a known integer number. We observe one path 
of the process at N discrete times, that (X(0),X(A N ),- ■ ■ ,X(NA N )). Let [f min , f ma x], with < f min < f max , 
be the chosen frequency band (see section 5, for an example). We discretize a (slightly modified) frequency band 
and compute the wavelet coefficients at the frequencies {fk)o<k<a N where 

fk = (qN) k for k = 0, • • • , a N , q N = ( an d a N = NA N . 

P \ fmin a J 

For notational convenience , we assume here that NAn is an integer number. By definition, we have /o = fmin/ 'P 
and f aN = fmax/ot, then, using the wavelet coefficients at the frequencies (fk)o<k<a N , we could detect all frequency 
changes (ui*) included in the band \f m im f ma x[- To simplify the notations, we use the following assumption : 

Assumption (C) : u* e]/ mm , f max [ for all i = 1, . . . ,K. 

In this framework, the estimation of the different parameters of X becomes a problem of linear regression with 
a known number of changes; thus, we follow the same method as in Bai (1994), Bai and Perron (1998), Laviellc 
(1999) or Lavielle and Moulines (2000) and define the estimated parameters (T^ N \ A^) as the couple of vectors 
which minimize the quadratic criterion : 

K+l tj + l-TN 

Q (N) (T, A) =J2 li'.- XiXj\ 2 ,and thus 

j=0 i=l+tj 

(fW.XW) = Argmin {q (n \T, A); T e A^ , A e B K } 

with 

• Y t = log(Jjv(l//i)), Xi = (log/,, 1) for i = 0,- • • ,a N ; 
log(/?/a)- 



• Tat 



log q N 



where [x] is the integer part of x. 



T = (i ,ti, • • • ,t K +i) G A [ p where 



A k -{(tor- ■ ,t K +i) G N K+2 ;t = 0,t K +i = a N + T N ,t J+1 -tj > t n for j = 0, • • • , K) ; 

/ — (2H- +1) \ 
A = (Ao, • • • , \k) G Br where Aj = I , 2 , i 3 v i i\ ) an d then 



log cr| + log K Hj (tp) 
B K ={{\ ,--- ,\ K ) with (H^a^eJC for all j e {0, 1, • • • , if }}. 

The integer t^t corresponds to the number of frequencies in the transition zones and log fi+ TN = log fi + log(/3/a). 
Obviously, for j = 0, • • • , K, the vector A^' provides the estimators Hj N ^ of H* and crj^ of a* by the relation 

/ -{2H (N) + l) \ 

= I , ( ,~.(m^n\ 3 , / n I ■ For a given T e .4^, each A^ is obtained from a linear regression 

1 y\og[(a) N) ) 2 }+\ogK^ N) {i,) J K 1 3 

of (Yi) onto (Xi) for i = tj + I, ■ ■ ■ , tj+i — tat. Thus, with T = {tj)a<j<K+i obtained from the minimization in T 
of Q( n ^(T,A), we define the different estimators of the change frequencies as 



..(AT) _ ^ ^ _ ^ /min f fmax P \ 

\ fmin & , 



^(iv) /. Jmin / Jmax H \ r ■ 1 7^ / nn \ 

= a /jtN) = a • — — ( — -) forj = !,••• , AT. (20) 



We have the following convergence : 
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be a discretized 



Proposition 4.1 Let X satisfy Assumptions (C) and (Bk) with a known K, (X& N , ■ ■ ■ ,Xna^ 
path, and ip satisfy Assumption (Al). Let Ajv be such as NAn — > oo and N(An) 2 — > when N — > oo. Assume 

, K . Then for all e > 0, there exists < C < oo such as for all large N , 



that (iif^of ) £ K. for all i = 



>C)<e for j = l, 



(21) 



Remark 4.2 The proof of this proposition shows a more general result, i.e. for (p, q) £ [3/4, 1] x [0, 1], for e > 0, 
there exists C > such as 



~(iV) 



■ K 



>Cj<e for j = l, 

with on = (NAjs[) q . For numerical considerations and convergence rate of the following estimators of the param- 
eters, we are going to fix now on p = 3/4 and q = 1 and then ajy = NAn- 

For j = 0, ••• ,K, the natural estimates of H* and a 2 * are given by the regression of (Yi) onto (log/j) for 
i £ $j N \ ■ ■ ■ , t^j+i ~tat}. But the probability that $j N \ tj+\ ~ tn] C [t*,t* +1 — tn] does not increase fast enough 
to 1 as N — > oo, in order to obtain a sufficiently fast convergence rate for these estimators. We address this 
difficulty as follows. We fix an integer number m > 3 and for j = 0, • • • , K, we consider [Uj N \ ] an interval 



strictly included in $j N — tn], such as 



L>i W) =e ) + 



t J+1 - tj - T N 



m + 1 



and vf°=*f> 



(22) 

Then we estimate the parameters from a regression onto m points uniformly distributed in [f7_- , vj ]; it provides 



AN) _ AN) 
m + 1 



TAT 



the following estimator \j from a regression of (Yi) onto (X,-) for 



i£{U) 



(N) 



V} N) }= W)"' + (k-l) 



(N) 



AN) AN) 

*J+1 _ *J - Tw 



m + 1 



= ( - (2jzf° + 1), log oar ' +log^ w (V) 



J J l<k<r. 
I 



(x^yxp 



and for all k = 1, 



5*« 



3 
Q 



aw 



, m, define g$(k) 

k/(m+l) 



fa 



(f)vvC) 



P \ fmin 

for all j G {1, — , if — 1} 




By this way, define 



(tog/* . 1 ) i6{ &w 



e{(7 3 (JV) ,-,v; <N) } 

f ^ fc/(m+l) 
Jmax 



frr 



and 



We get the following central limit theorems for the corresponding estimators , a 2 ^ "* ) 

Proposition 4.2 Under the same assumptions as in Proposition^. 1\ for all j = 0, • • ■ , K , 



{NA N f' 2 ( 



(JV) 



x> ^ 



AT(0,r^) 



w/iere I\ 3 = ( X* X* 
lowing matrix : 



X*Y,*X* X* X* 



with X* = (log ff ; (A) , i) 1 



<fc<r 



(23) 



( s £-z)l<fc,Z<™ i/ie 



*j 

S kl — 



2H' 



2 ■ 




1? 



(24) 



Remark 4.3 Another possible choice would be to consider the regression for all the available frequencies in the 
interval [Uj , Vj ]. The number of considered frequencies increases then with the rate on — NAn- However, it 
does not improve significantly the convergence since the remainders of the regression are very strongly dependent. 
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4.2 Goodness of fit test 

It is also possible to estimate parameters H* and tr* from an feasible (or estimated) generalized least squares 
estimation (for more details, see Amemiya, chap. 6.3, 1985). Indeed, we can identify the asymptotic covariance 
matrix E* for j = 0, • • • , K : this matrix has the form E* = E(iT*, oj*, and, from the previous limit theorems, 

E^ = T,(Hj ,Wj N \u}j^l) converges in probability to E*. Thus, it is possible to construct an estimator Xj N ^ of 
Xj with a feasible generalized least squares (F.G.L.S.) regression i.e. by minimizing 

|| y} n) - xf h ||| ( „,= (y w - xf hy (sf >) 1 (ff > - Xf >A). 



First, we give asymptotic behavior of A,- 



-(2H^> + l),loga^ +log^,(V) 



3 



Proposition 4.3 Under the same assumptions as in Proposition for all j = 0, • • ■ , K , 

(NA N )V 2 U N ) - X*) -H* AT(0,r?) (25) 



with T 2 " 



: -v;) : . 



For ? = 0, • • • , X, the vectors Y} N ^ and xi^A^ are two different estimators of the vector 

(-(2H* + ljlog/i + log of + log^.(V)) ie rr/(w) !>(«)!■ It suggests to define the following goodness of fit test, 
is defined as the su 

frequency ranges: 



The test statistic is defined as the sum of the squared distances between these two estimators for all K + 1 




if = (NA N ) • | £ || yW - xf)Af ) 



This distance is the F.G.L.S. distance between points (log/*, ^) ie {;y( N ) y (JV) } f° r J = 0, • • • , A an d the (if + 1) 
F.G.L.S. regression lines. As a consequence, we get 

Proposition 4.4 Under assumptions of Proposition \j~T\ we have 

T [ K N) X 2 ((A+l)(™-2)). (26) 

Remark 4.4 Proposition \4-4\ niay be explained with heuristic arguments. Remainders are turned white, thus it 
is only natural for the sum of the second regression remainder squares to asymptotically form a \ 2 process. The 
number of degrees of freedom is (K + l)(m — 2) because one loses two degrees of freedom after the twice estimation 
of the (K + 1) vectors A* (we also show that these vectors are asymptotically independent). 



4.3 Estimation of the number of frequency changes 

Throughout the previous study, the number of frequency change, A, is assumed to be known. But the previous test 
provides a way for estimating A. In fact, it can be recursively done by beginning with A — and continuing till 
the assumption "A is a (M^)-F.B.M." is accepted. The following applications in biomechanics provide different 
examples of the power of discrimination of such a procedure. However, this estimation of the number of frequency 
changes must be carefully applied : from numerical and heuristic arguments, it does not seem reasonable to work 
with A > 2. 
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4.4 Estimation procedure and on the choice of parameters 

Thus, for identifying a (Mr )-multiscalc fractional (with K unknown) from a time series (Xo, Xa n , • • • , Xna n ) we 
suggest the following procedure: 

1. Begin with K = 0. 

2. Choose a mother wavelet -0 (and thus a and /3), a frequency band [f mini fmax] and to (see below for these 
different choices). 

3. Compute the different frequencies (fi)o<i<a N ■ 

4. Compute the vector (Y l ) < i < aN = (log JnO-/ fi))o<i<a N ■ 

5. Minimize Q^ N \T, A) and thus compute the different values of wj JV ' 1 for j = 1, • • • , X. 

6. Compute the different regression moments , • • • , Vj } and then the estimators Xj (for j = 0, • • • , K). 

7. Compute the different matrices £j ^ and then A~- (for j = 0, • • • , if). 

8. Compute and compare its value to the 95%-quantile of a x 2 ((AT + l)(m — 2)). If the test is rejected 
then go back to step 2. with K = K + 1. 

How to chose the function ip and the parameters f m in , fmax and m ? 

1. Choice of ip : The mother wavelet i/j has to satisfy Assumptions (Al) but as we say previously it is not 
mandatory to associate this function to orthogonality properties. However, the Lemarie-Meyer wavelet is 
a natural choice with good numerical properties of asymptotic decreasing but a too large ratio (3/a which 
implies a too large transition zone of frequencies. The function ip can also be deduced from an arbitrary 

construction of its Fourier transform ip; for instance, we propose ipi(X) = exp (77—; r--z r— p- ] l a <\x\<3 

\{\X\-a)((3- \X\)J -' 

and the function built from a translation of the Fourier transform of the Lemarie-Meyer function to 
[— 2ir, —ir) U [k, 2tt] (thus the ratio is now (3/a = 2). The results obtained from those functions ipi and tp2 are 
essentially the same than with the Lemarie-Meyer mother function, they appear more precise for the detection 
of frequency changes u* (because log/9/a and thus the transition band, could be as small as wanted) and 
less precise for the estimation of parameters H* (because ipi and ^2 are not concentrated as well around 0). 

2. Choice of f m %n an( ^ fmax '• (we assume here that the frequencies are given in the inverse of (X±, X2 • • • ) 
time unity). The choice of f m in and fmax is first driven by the selection of a frequency band inside which the 

process has to be studied; the inspected frequency band is then [ ™" , mw l. Secondly, N x ™" should be 

pa p 

P fmin 

large enough for computing In(~z ) m ©■ Formally one only needs to have N x — - — > 1 but numerically 

Jmin P 

N x > 10 seems to be necessary to use correctly the central limit theorem. Finally, the discretiza- 

P 

(X 

tion problem implies that fmax cannot be too large for providing a good estimation of dx(— f j^Aat) by 

J max 

ex (7 ,kAx). In practice ^ max < _ — appears as a minimal condition. 

fmax OL Ajv 

3. Choice of m : Formally, m could be chosen such as 3 < m < minj(t* +1 —rjf — tj). Theoretically, the larger 
the to, the closer to 1 the power of the test. But numerical considerations imply that if to is too large then 
the different matrix £j- are extremely correlated and the quality of the test is very dependent to the quality 
of the different estimations of A*. As a consequence, we chose 5 < m < 10. 
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5 Numerical simulation and applications in Biomechanics 
5.1 Simulations 

1. Initially, we apply the estimators and tests to several simulated trajectories of classical F.B.M. (generated 
according to a Choleski decomposition) with different values of H = 0.2, 0.4, 0.6 and 0.8). The selected of values 
of different parameters are : N — 6000, = 0.03, m = 5, / m j n = 0.05 and f ma x = 20. There are 30 independent 
replications of each time series. The results are presented in the following table : 



Theoretical values of H 


0.2 


0.4 


0.6 


0.8 


Empirical mean of H 


0.148 


0.384 


0.599 


0.821 


Standard deviation of H 


0.034 


0.031 


0.041 


0.048 



The Figure 4 presents the log-log representation for one trajectory : the linearity is seeming. Moreover, Figure 5 
exhibits a histogram of the distribution of the test statistic T(j (in this case K = and 30 x 4 = 120 independent 
realizations) compared to a x 2 -distribution with 3 degrees of freedom. The goodness-of-fit Kolmogorov-Smirnov 
test for to the x 2 (3) distribution is also accepted (with D ~ 0.091 and p — value ~ 0.272). 




Figure 4 : The log-log representation for a trajectory of a (Mo)-FBM (left, with H — 0.6) and (Mi)-FBM (right) 

2. Then, we apply to 30 independent replications trajectories of (Mi)-FBM (generated according to a Choleski 
decomposition with numerical approximations of the covarianccs ) with H = 0.2 and <jg = 10, Hi = 0.7, and 
of = 5, and U\ — 5. The results (with parameters : N — 6000, = 0.03, m — 5, / m j„ = 0.8 and f max = 16) are 
the following : 



Theoretical value 


H = 0.2 


Hi = 0.7 


u>i = 5 


Empirical mean 


0.197 


0.693 


5.18 


Standard deviation 


0.110 


0.068 


0.491 



Figure 4 presents the log-log representation for one trajectory, with the 2 regression lines. The hypothesis of the 
modelling with a simple FBM (therefore with K = 0) is always rejected (in such a case, the model is misspecified 
and the statistic Tq N ^ is then between 39.3 and 126.8, very different from the realizations of x 2 -distribution with 
3 degrees of freedom). On the contrary, the hypothesis of the modelling with a (Mi)-FBM is always accepted and 
a histogram of the realizations of the test statistic is presented in Figure 5 (compared to a x 2 -distribution 

with 6 degrees of freedom). The goodness-of-fit Kolmogorov-Smirnov test for T x to the x 2 (6) distribution is also 
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Figure 5 : The empirical distribution of T and T{ (respectively) compared to the corresponding % 2 distri- 
bution in the cases of simulated trajectories of (Mo)-FBM (left) and (Mi)-FBM (right) 

Conclusion of these simulations : the results are surprisingly good compared with the complexity of the 
method. The asymptotic distribution of the test statistics can be used for real data. However, the computation 
time is important (especially for the computation of the test statistic) : 3 hours are necessary for the treatment of 
each (Mi)-FBM replication. 

5.2 Applications in Biomechanics 

We apply our statistics to different trajectories (see the description in the Introduction) with the following param- 
eters : 

• N = 6000 and A N = 0.03; 

• The mother wavelet is il>i (with a = 5 and f3 = 10). 

• The choice of the frequency band is f m in = 0.15 and fmax = 15 which corresponds to a detection frequency 
band [0.52 , 38.32] Hz (with mother wavelet tpx); 

• m = 5. 

First, we study the X-trajectories of one subject (fore-aft direction) for different feet position (0; 2; 10; 20 cm 
clearance and 0; 15; 30; 45° angle). In all the cases, the test (with a type I error of 5%) rejects the hypothesis of 
a modelling with a simple (Mq)-FBM. But the modelling with a (Mi)-FBM is accepted by the test 12 times out 
of 16, with an empirical mean of lj\ ~ 3.5 and a standard deviation of u5i ~ 1 (the different values of Hq and H\ 
are in [0.9, 1] in the different cases). 

For the different y-trajectories (medio-lateral direction) of the same patient, the test rejects the hypothesis of 
a modelling with a simple (Mo)-FBM in all the case. The modelling with a (Mi)-FBM is accepted by the test 13 
times out of 16, with an empirical mean of cDi ~ 3.1 and a standard deviation of lj\ ~ 1 (the different values of Hq 
and Hi are in [0.8, 1] in the different cases). 



Figure 6 presents log-log plots of Jjv(/fe) versus fk (i.e. log Jjv(/fc) vs. log fk) of all the experiments, for X- 
trajectories (left) and Y-trajectories (right). 
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Conclusion of these applications to biomechanics data : all these results allow us to give new interpretations 
on the upright position. The behavior of X-trajectories and Y trajectories are very similar, for all the positions 
of the feet (the studied statistics do not seem to depend on the angle and clearance of the feet). The (Mi)-FBM 
models of these trajectories fit well, which suggests two different type of behavior for low and high frequencies. 
The frequency change is around 3 Hz, which corresponds to a physiological change : this could be interpreted 
for instance as the passage of a cerebral control of the stability by the inner ear to a muscular auto-stabilization. 
We return to [12] for a more detailed discussion of the biomechanical interpretations. Such an estimation of this 
frequency change would be very interesting for a better detection of certain pathologies and to help in their cure. 

A Proofs 

A.l Proof of Theorem 13.11 

First, we prove the following technical Lemma : 

Lemma A.l Let X be a (M K )-MBM. For (t,f) e R 2 , and (u,u') e R\, define : 

S(t, u, t', v!) = E \(X(t + u)- X{t)) ■ {X{t' + v!) - X . (27) 

1. For all (u,u' ,t,t') G IR 2 ^ x IR 2 , there exists a constant C > depending only on the parameters (uJj)j, (o~j)j 
and (Hj)j such that : 

\S(t,u,t',u')\ < C ■ (u Hk ■ l u <i + u Ho ■ l u> i) x (u> Hk - l u ,<i+w Bo - V>i); (28) 

2. More precisely, if (max(«,u') • ojk) < 1 o,nd max(n, «') < — • \t' — 1| , there exists a constant C > 
depending only on the parameters {oJj)j, and (Hj)j such that : 



\S(t, 2u, t>, 2u')\ <C-(u.u> + max( U , .') 4 ) + M , _ M | + { \t -t> + J - u\>->* }) 

Proof. 1/ First, the Cauchy-Schwarz inequality implies that 

\S(t,u,t',u')\ < yjE[(X(t + a)- X(t)Y] x ^E[(X(t' + a')-X(f) 2 )]. 

2 J^-r f auJ j+l fl — COSW) 

But, E [(X(t + a) - X(t))] =4 22 a 2 ■ a 2H > / v 2g +1 ; dv. Then, the following expansions : 

j=0 J au) j v 



(29) 



x (i_cos^) dv = I ^^)^ 2H + °(^ 2ff ) fOT 
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with C(H) 



(1 — cosu) 

— czu, imply that : 



£[(*(t + M )-x(t))] : 



that achieves the proof of the majoration 



4 • o% ■ C{H K ) ■ u 2H '< + 0{u 2 ) when u -> 0; 
4 • erg • C(.ffo) • w 2ffo + 0(1) when u -> oo; 



(30) 



2/ We turn now to the proof of the upper bound (|2"5|) . To begin with, we remark that for all (t, t' , u, u') S 5? 4 , the 
following equalities are true : 



S(t,2u,t',2u') = 



-i(t+2«)£ 



t£\/ e i{t'+2u')£ _ p it'C 



P 2 (0 



d( 



p 2 (0 

sin(u0 • sin(u'0 • cos (£(t' — t + v! — u)) 



p 2 (0 



d£ 



N - . , sin«) ■sin(-u^) ■ cos (£(f - f + u' - u)) 
Z^ CT * / ■ ?2iS+i 



i=0 



(31) 



Then, we bound the different integrals. 

• First, we threat the case i = K that is when the upper limit of the integral is oo. In this case, we can rewrite 
the integral between ojk and oo as the difference of the integral between and oo and the one between and ujk, 

poo poo r^K 

that is / 7(0 <£= / 7(0 d£ - 7(0 df where 7 (0 = sin(<) • sin(z/0 • cos (£(i' - t W - a)) r (2HR " +1) - 

Jujk Jo Jo 

The second integral of the right hand side can be bounded by the same argument than the terms of order i = in 

(f3Tj) . The first one corresponds to the expression of the covariance of the increments of a F.B.M. Bh k with Hurst 

parameter Hk and variance 1. Thus, for all (t, t') € M 2 , (u,u r ) € JR^_ such that 4 • max(ti,«') < \t' — t\, we get : 



sin(u0 • sin(u'0 • cos (£(t' — t + u' — u)) 



E 



2C 2 {H K ) 



£2H K + 1 

({B Hk (t + 2«) - B Hk (t)) • (S ffK (t' + 2u0 - B Hk (f ))) 

i - t' + 2u| 2Ha ' - It - t' + 2n - 2u'\ 2Hk - It - t'l 2 ^ + It' - t - 2u'| 2Hk 



1 



|t - t' + u' - u\ 



1 -j/ 1 2 — 2Hk ' 



with D(Hk) > 0. 



(32) 



• Next, we consider the integrals with a finite upper limit and a non-zero lower limit. This corresponds to 
i = 1, . . . , K — 1. In these cases, for b > 0, an integration by parts provides us 



m+1 sin«) • sin(u'0 • cos (60 _ 1 

flc; t" 



i2Hi + l 



u • cos(iif) • sin(u'0 + u' ■ cos(w'0 • sin(w0 ,, iS 

^-^ • sm(60 d£ 



£2H z + l 

sin«) • sin(u'0 • sin(60 1 



2Hi + l 



(33) 



By using the majoration |sin(ua;)| < ux, |sin(w'a;)| < u'x, |cos(itx)| < 1, |cos(u'x)| < 1 and |sin(6a;)| < 1 for 



x > 0, we deduce that for all (u, u', 6) G 

ru,i+1 sin«) • sin(u'0 • cos(60 



2Hi + l 



< Ci- 



u • u 



b ' 



(34) 
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where Q > is a constant depending only on Hi,Ui and Wj+i. 

• Finally, it remains to bound the two integrals with lower limit 0. We will show only how to bound 

r 1 sin«) • sinfu'O ■ cos(6£) , , . , ml . . , , 

/ g a£, since the other integral can be treated similarly. The integration by part formula 

(f3"3"| remains valid even when the lower limit is 0. Indeed, the integrand can be bounded by C x ^^ 2Hi and 

^i-2H< ^ <; oo ^ soon as Hi < 1. After this remark, we bound the three terms of the right hand side of (|33|) . 

i) From | sin(itx)| < ux, \ sm(u'x)\ < u'x and | sin(te)| < 1 for x > 0, we deduce that for all (it, v! , b) £ M+, 

sin«) -sm(yO -sin^) ]" 1 , ,i-2H \ „ j 

S Jo 

ii) For all (£,£') £ [0, the power series expansion of x i— > sin(ai) implies that 

00 fc ,23+1 . ( u A2(fc-j)+l 



(35) 



sin«) • sin«) = u • u' • C 2 + J! ( £ ( 2 "+ i 



fc=l 3=0 



(2j + l)!-(2(fc-j) + l) 



im)-^ 



fcf2fe+2 



One can remark that 

fc 

E 



^ o (2i + l)!.(2(fc-j) + l)! " max ^ u > Z 5(2i + i)!.(2(fc_j) + l)! 



< max(it, u ) 



/\2fc+2 



because X^>o (2j+"IJi — ^- As a consequence, when (max(«, 11') ■ uoi) < 1 and b > 0, integration and summation 
can be interchanged and 



Jl sin«) • sin(u'£) • sin(&0 



2H0+2 



d£ — u ■ v! 



sin(k£) 



\2k+2 



00 / 

— ( maX (' U ' 

fc=l ^ 

/\4 , ,3-2H 



2fc-2H 



But 



sin(bQ 
l 57 ^ 



2-ffn-l 







< 



b ' Wl sinffl 

d£. Denote M(ff) = sup 



max(u, it') ■ u>i~ 
1 — (max(u, u 1 ) ■ cji) 2 

sin(0 



xeM 



sin(u£) • sin(u'£) ■ sin(6£) 

£2H +2 

iii) Similarly for (£,£') S [0,wi], we have 



£2H 
3-2H, 



< M(H { 



( ^ 1 + ^n^ 
1 — (max(zt, u'J • 



for < H < 1. Thus 



(36) 



00 re 

cos«) • Bin(u'0 = u' • i + ( £ 



fc=l j=0 



u 2j . („/j2(fc-i)+l 

(2j)!-(2(A:-j) + l)! 



(-ire 



fc e 2fc+l 



But for k > 1. 



re 

3=0 



u 23 . ( u /)2(k-i)+l 



(2j)!-(2(fc-j) + l) 



— ^ < max(u, w') 2,£+1 . As a consequence, when (max(u, it') • u>\) < 1 and 



& > 0, integration and summation can be interchanged and we get 



u ■ cos(w£) • sm(u'£) • sin(&£) 



2H +i 



d^ — u ■ u 



sin (b£) 

£2H„ 



d£ 



r\4 , ,3-2H 



< 



max(u, u') ■ uj 1 
1 — (max(u, u') ■ loi) 2 



(37) 



Therefore from |35|) , (|36|) , (|3T[) and (|32| . we deduce for (u,u') such that max(«, «') ■ < 1/2 



sin(u£) • sin(u'^) • cos(6£) 

£2H K + 1 



< (£>(ffjf ) + 3M(H K )) 



b 2 



-2H h 



4- 



max( M , mQ 4 3 _ 2ffK 



1-2H A 
K 



By combining the two previous bounds with (|3 1 [) and (fM)) . we deduce (|29p and this finishes the proof. 



The proof of Theorem 13.11 uses the two following lemmas: 
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Lemma A. 2 Under the same notations and assumptions as in Theorem \3.1[ there exists two constants C\ > 
and G<i > depending only on r, a m i n and a max such that for all N 

i) sup max Ee 2 N (a,k) < d x ip(N) (38) 

ae[a m in,a max ] k£D N (a) 

ii) sup max ly/ai £Ar(ai, fc) — v / a^e A r(a 2 , fc)| 2 < C*2 x (p(N) x |a 2 — ai | 2 (39) 

oi,a2e[o m i„,o mM ] fee-DN(ai)nDjv(a 2 ) 

w/iere p(iV) = + A]+ 2Hk + A 2 N log TV + Aj+^AT"^^ + (iVA^)" 2 with H = max{H, , i = 0, • • • ,K}. 
Hi) Moreover (N Ajv) <p(N) — > w/ien iV — > oo. 



Proof. The error £jv(a, fc) contains three different terms, the first one corresponds to the replacement of the 
integral onto the interval [0, T/v] by its Riemann sum, the second and the third ones correspond to the replacement 
of the integral onto Ft by the integral onto the interval [0,T/v] where Tn = NAn- More precisely, we have 

e N {a,k) = — x (e 1N (a,k) +s 2 ,N(a,k) +£3,N(a,k)) (40) 



with 



e hN (a,k) = / ^(--kA N )X(t)dt - A N J2^(^ JL -kA N )X(pA N ), 
Jo ^ 



f°° t 

£2,jv(a,fc) = / ip(- -kA N )X(t)dt, 
Jt n a 

f° t 

£ 3 ,jv(a,fc) = / Tjj(--kA N )X(t)dt. 

J-oo a 

By using (x + y + z) 2 < 3 (x 2 + y 2 + z 2 ) for all real numbers a;, y, z, we deduce 

/ 3 \ 



Es 2 N (a,k) < (jj) x [Y,Esl N {a,k)\ . (41) 



=i 



We now bound the different terms E e 2 N (a, k) for i = 1,2,3 : 
(1) Bound of Eel N (a,k). 

We have the decomposition £uv(a, /c) = /^(a, k) + /2.w(a, k) , where 

tf-i o( P +i)A N t 



h, N {a,k) = V/ ^(--fcAjv) [X(t)-X(pA N ))dt 

and h,N(a,k) = J2 J (1>{~ ~ kA N J - - kA N j J X(pA N ) dt. 



P =o j p a k 

Then, the inequality (x + y) 2 < 2(x 2 + y 2 ) which is valid for all (x,y) G M 2 , implies 

E el N (a, k)<2E (I* N (a,k))+2E {l% (a, k)) . 

On one hand, we have 

AT-1JV-1 r (p+l)A N ,(p' + l)A N t t , 



E(lf N (a,k)) = Y^J2 / dtdt , V(--*A J v)v(--A;A J v)£7 

w-uv-i ,a«/-a„ ,., , .„ x ,,,, 



(X(t) - X( P A N ))(X(t') - X(p'A N )) 



EE/ 7 dudu'jj^^ - kA N )j;( U ' + P ' AN - kA N )s{pA N ,u,p'A N> u'). 
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Afterwards, we use Lemma IA.1I to bound the terms S(pA^, u,p'A^, u'), where we use different types of bounds 
depending wether is in the vicinity of the diagonal or not. Namely, when \p — p'\ < 3 we use the upper 

bound l|28p , and otherwise we use the one in (f2"5|) . Observe that the assumptions of Lemma IA.1I are satisfied for 
large enough N since NAn — > oo, as N — > oo. Thus, when N is large enough, we get 



E(l 2 N (a,k)) < J J dudu' 4>{ 



p = O.p' = 

Ip-p'I < 3 



' u + pA 



N 



kA 



A 



u' +p'A N 



kA 



A 



C-{u- u') Hk + 



N-l ^AjyAjv 

/ / dudu 1 

„ , Jfi Jo 

p = 0, p =0 

\p — p > 4 



u + pA 



N 



kA N )^ U ' + ^ AN -kA N 



■■ x C ■ A 



s 1 (\p'- P \-i)A N 1 i^r,K\((\p'-p\-i)A N ) 2 - 2 ^ 



N-i / a i \ 

< CAl +2H « £ SU P ^ (— -W 



p = 0,p' = 

Ip-p'I < 3 

JV-l 



x sup 

e'e(o,i) 



+C-A^ V sup VI ( 

ee(o,i) \ v 



p = 0,p' = o 
|p-p'| > 4 



k)A 



N 



x sup 
0'e(o,i) 



fc A 



A 



fc A 



1 



max 



(|p' - p\ - 1)A N i=o,- ,x I ((b' - p| - 1)A N ) 



1 



\2-2H; 



However, according to Assumption (Al), for every integer m E IN* there exists a constant C > such that for all 
x e M, \ip(x)\ < C ■ (1 + |x|)~ m . In particular, sup 

0'£(O,1) 



is bounded. Therefore, 



N-l 



E(l 2 N (a,k)) < C-A 2 + 2Hk Y\ S UP (l + \0+p-ak\A N /a 

p=o <>e(o,i) v 



iV-l 



A" 



=0 «e(o,i) 



C-A^^ sup (l + |0+p-afc|Ajv/a) ^ 



q=3 



qA N i=o,-,K { (qA N ) 



1 



)2—2Hi 



But 



jV-1 



A N V sup (l + lfl + p-ofcl — ) < CA,v V (l 
p^ 06(oa) V a ; ^ 1 



_- — OO 

OO 



A 



< c 



1 + — ) dx 

a 



= C\a\ {l + \y\)- m dx < C\a max \. 



Let us denote H = maxjff, , i = 0, • • • , X} and H_ — min{7Ji , i = 0, • • • , i^}. We deduce 



S (7 1 V(a I fc)) < C-A^ + C.A^f:(^ + i ma, K {^^}) 

g— 3 



C A 



1+2H K 
N 



A 2 N l giV + Ai r (-^XJ- 3 i 

< C (A™ + A 2 ^ logiV + A 2 N + A^N-^ 2 " 



1 * 1 

— y — 

2-2H l—< „2- 



N g=3 y A q=A~ 1 ^ 



2-2H 



(42) 
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On the other hand, by using Lemma I A. 11 formula (|27p . we get 

N-1N-1 r (p+l)A Nr (p'+l)A N / 

£(lf >iV (a,fc)) = EE / / dud«'U(--*A w )-V( 



AT 



fcA 



A 



^-fcAjv) -</>( 



y ' Ax fcAiv) ) 5(0,^,0, p'Ajv) 



' u + pA 



< C 2_, / du ■ — ■ sup 

(AT-1 
V sup V'( 
p=0 te[o,A„] 



A' 



fcA 



TV 



ft 

a 



fcA 



a 



x |(pAjv) ffo + ( P A N ) 



H K \ 



,/t + pA N 



<//( 



fcA 



x \(pA N ) H " + ( P A N ) 



H K \ 



,/t + pA 



A 



fcA 



A 



x \{ P A N ) Ho + { P A N ) 



H K 



But Assumption (Al) implies that for m = 4, there exists a constant C > such that for all x E JR, \ip'(x)\ < 
C-(l + |x|)- m . We deduce 



N-l 



p=0 



,/t + pAi 



Aat > Sup 1p I fcAjv 

- te[o,A^i v a 



<CA N E 



x |(pAat) h ° + (pA N ) HK 
1 



(l + bAjvl 



00 \x\ H ° + |a;| g « 
-co (1 + M) m 



dx < 



Therefore, 



E(ll N (a,k)) < CA%. 



(43) 



(2) Bound of Ee\ N {a,k). 

By using Lemma I A. 1 1 and Cauchy-Schwartz inequality, we deduce that for N large enough, 

/>oo />oo / / \ 

£^ >JV (a,fc)) = y y V (- - ^AjvJ ^> ^— - kA N j S(0,u,0,u')dudu' 



IT N 

' i N 



- fcA 



N 



,2H, 



°du 



On one hand, fc € Dn(cl) implies that k < [(1 — r)iV/o]. On the other hand, u > Tjv = A" • Ajv- Therefore, we 



have ( 1 



(l + \^-kA N \) 



> (u — (1 — r)iVAjv) /a. This implies that for m > 4 and N large enough, 

\ — m 

to 

dv, 



fcA 



A 



— m 

,r""d„ • / (--(l-r)^M -/-'" >du 
It n \a a J 

r ~ (« + 1 - r) 2ff « 



ft" 



(JVA 



im-2ifo-l 



U" 



by making the change of variable u — (NAn) (v + 1 — r). Consequently, 

Eel jN [a,k) <C ■ ' 



(iVAjv) 2 



(44) 



(3) Bound of ^£3^(0,^;). 
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kA N 


<-< 


/> 






a 



By using the same kind of argument than in (2), one obtains that for N large enough 

Eel N {a,k)) = J J t/> (~ - fcAjy) tp f— - kA^ S(0, u, 0, v!) du du' 

)—m 
(\u\ 2Hk + \u\ 2H «)du 

< C (r ^(v~r/2.NA N?d Y 

\Jr/2-NA N vm J 

As a consequence, for m > 4 and N large enough, 

(45) 

Finally, from (O, (O, (gU) and (gHJ), we deduce that ([35]) holds. This finishes the proof of the point i). Since 
NAn — * oo and iVA^ — * 0, (NA]y)(p(N) converges to when N — > oo. This proves the point iii). To complete 
the proof of Lemma [A.2l it remains to proves the point ii). We deduce from the decomposition (|40[) that 

3 

E |-v/aisjv(ai) k) - % /a 2 "£Ar(a 2 , k)\ 2 < 3 E \e^ N {a 2 , k) - e i>N (a 1 , k)\ 2 (46) 

i=l 

The same calculations than the ones used to prove the point i) provide the upper bound on the terms 

E |£i,jv(a2, k) — £ijv(ai, k)\ 2 . Indeed, consider for instance the terms with £2, at, then by using Taylor formula, for 

every pair (ai, a 2 ) with a m i n < a\ < a 2 < a max there exists a real number 9 € (ai, 02) such that 

£2,jv(a2,fe) -e 2 ,N{ai,k) = (a 2 - 01) x ^ ^-^^ -0' Q - fcA^J * 

Next, by using the same kind of arguments than for the bound of i?£ 2 jv( a j k) in point z,), we get that for every 
integer m > 4, every 01, <z 2 in [a mm , a ma J and S D N (ai) n D N (a 2 ) 

£ |£ 2 ,jv(a 2 , fc)-£ 2 ,iv(ai, fc)| 2 < -^|a 2 -ai| 2 / (1 + ~ - fcAjv ) U 1+2Ho du 

< C\a 2 - ai \ 2 (NA N )- 2 . 

We deduce similarly that 

E \e 3 M(a 2 ,k)~e 3 M(a 1 ,k)\ 2 < C \a 2 - ai \ 2 {NA N y 2 . 
At this point, it remains to show 

E \e hN (a 2 ,k)-e hN ( ai ,k)\ 2 < C \a 2 - ai \ 2 <p(N) (47) 
to finish the proof of item ii). But, we have the decomposition 

E |ei,jv(02, k) - £i,jv(oi, k)\ 2 < 2E |/i,jv(o2, k) - h,N(ai, k)\ 2 + 2E \I 2 , N (a 2 ,k) - 7 2 ,jv(ai, k)\ 2 
However, Taylor Formula implies the existence of two real numbers 9\, 9 2 £ (a%, a 2 ) such that 
h,N(a 2 ,k) - 7 1)JV (ai, k) = (a 2 - a%) ■ Ii, N (9i,k) for i = 1 or 2 

where Ii t pf(a,k) is obtained by replacing into the expression of Ii^(a,k) the map -0^ kAj^j by the map 

Pj X^- fcA w ) and - fcA^v) by x - fcA^v). So, 

-B |ei,jv(o2, fc) — £i,w(ai, fc)| 2 < C |a 2 - ai| 2 x {£ 7^(01, fc) + E% N (0 2 ,k)} 
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Since the map t i— > f — j x ip' y fcAArJ is still continuously differentiable and fast decreasing, one can lead same 

calculations that in the bound of EI 2 N (a, k). We finally get E if N (02, k) + E I^ N ((>2, k) < Ctp(N). This implies 
(|47f and completes the proof of Lemma [A.2l ■ 



Lemma A. 3 Under the same assumptions as in Theorem \3.1\ there exists a positive constant C > such that for 
every real number a > and N € IN*, we have E |//v(a) — J/v(a)| < C x ip(N) 1 ^ 2 . 

Proof. Since the variables d = d(a, fcAjv) and e = e(a,kAN) are Gaussian, the variables d 2 — e 2 have finite 
second order moment and Jensen's inequality implies 



E \I N (a) - J N (a)\ = E 



]T (d 2 (a,kA N )-e 2 (a,kA N )) 

1 N(a> ' keD N (a) 

^ TrTT~\\ E \I E (<P(a,kA N )-e*(a,kA N )) 2 

Then we derive an upper bound for the expectations E (d 2 (a, A;Ajv) — e 2 (a, /cAjv)) 2 . Indeed, d and e are jointly 
Gaussian variables with zero means. One has 

E (d 2 - e 2 ) 2 =E(d- e) 2 (d + e) 2 =: E e 2 ' Z 2 , 

where e = d — e and Z = d + e are also jointly Gaussian and have mean zero. By using that Z = a<x er-f 1 pe + £, 
where er 2 = E e 2 , c 2 = E Z 2 , p — corr(s, Z) and where £ is independent of e and Gaussian, one can show that 

Ee 2 Z 2 = (Es 4 )^- + Ee 2 Ei 2 = Sa^p 2 + a\o\(\ - p 2 ) < 3afa 2 2 
a i 

= 3Ee 2 x E{d + e) 2 . 
But (d + e) 2 = (2d - e) 2 < 8d 2 + 2e 2 , therefore 

V6 



E\I N (a)-J N (a)\ < V £ ^Es 2 N (a,k)xy/E [id 2 (a,kA N ) + e 2 N (a,k)] 

< E x ' E E [4d 2 (a,kA N )+e 2 N (a,k)] 



keD N (a) ) {kED N (a) 

1/2 , >. 1/2 



< 



where the two last inequalities follow from Cauchy-Schwartz inequality and Ed 2 (a, kAjy) = T\[a) for every integer 
k. Since ip is compactly supported, then sup |2i(a)| < oo. By combining this remark with Lemma lA. 21 i). 

this provides E \In(o>) — Jn(o>)\ < C x ip(N) 1 ^ 2 and finishes the proof of the lemma. ■ 
Now, the following proof of Theorem 13.11 can be established : 

Proof. [Theorem 13.1] From Lemma [A.3I combined with Lemma [5T2] iii) , we deduce 

lim [N A N ) 1/2 E \I N (a) - J N (a)\ = 0. (48) 

N— ►oo 
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Combined with @, this implies the convergence of the finite-dimensional distribution in (|16p . Indeed, it suffices 
to show that 



y/NA N (log Ma) ~ log/jvO*)) 0. (49) 

Let e > 0. By using the inequality |log(x) — \og(y)\ < 2 \x/y — 1|, valid for all |x/y — 1| < 1/2, x, y > one can 
show that 



P I v N A N \\ogJ N {a)-\ogI N (a)\>e 



- ^ ! 2 \M«)/M*) - II > 7^^) +P (l-W/W - 1| > ^ 
<2P(|J A r( a )-/ 7V ( a )|>/^L) (50) 



2 V^VA 



<2Pl \ J N (a) - I N (a)\ < £ Xl( ° ) and | Jy(a) - J^(o)| > £/jv ^ 



eli (a) 



+ 2Fp)-/,(a)|> 4VM _ 
< 2F (|Jiv(a) - I*(a)| > ^==) +2P (l N (a) < ^) . (51) 

The second inequality in (f5D|) is valid for all N such that e/^NA^ < 1/2, that is, for all sufficiently large N. The 
second term in the right-hand side of (fSTj) vanishes, as iV — > 00, because /at (a) — ► 2i(a). By using the Markov 

N— >oo 



inequality, one can bound above the first term in the right-hand side of (|5ip by 



In [8], one have proved the tightness and the weak convergence of 



8y/NA^ 

£ JAr(a - /jv fl • 

eli (a) 

Thus, from (|48|) , one obtains Relation (|49|) . which completes the proof of the convergence of the finite distributions. 
To finish with the proof of Theorem 13.11 we have to show the tightness of the sequence (Ljv(a)) a < a < a 
where Lpf(a) — y/ N Ajv {Jn{o) — liip)^. Observe one has the decomposition Ljsr(a) = Lijv(a) + L2,Ar( a ) with 
L hN (a) = ^NA^flNia)-!!^ 

L 2 , N (a) = VWAH U N (a) - I N {a)^j 
(Li i Ar(a)) a < a < a m Skorokhod topology on the space of cad-lag functions on [a m i n , a max \. From the other 
hand, (|48p implies that L 2 ^{a) — ► for all a £ [dmin-j a>max\' Note that the limit process is null, thus it is 

N^oo 

obviously continuous. Then, provided one have shown the tightness of (L 2j i\r(a)) a < a < a 1 on e can deduce the 
tightness of (LN(a>)) a <a<a > see f° r instance Jacod and Shyriaev, Cor 3.33, p. 317. Next, one deduce the 
weak convergence of L^{a) to Z(a) in the Skorokod topology on the space of cad-lag functions on [a m i n ,a max ]. 
The last step is the proof of the tightness of {L 2 ^(a)) a . <a<a ■ Following Ikeda and Watanabe, Th.4.3 p. 18, 
it suffices to show the existence of a positive constant M 2 such that for all 01,026 [a m i n , a max ] 



E {L 2tN {a 2 ) - L 2>N ( ai )) < M 2 \a 2 ~ ai \ 2 . (52) 



Nowever, from ([5]) and (fF4"|) , we get 



J N (a) - I N (a) = \D N (a)\ 1 £ (e 2 (a, k) — d 2 (a, k)) 

keD N (a) 
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Therefore, for a\ < a 2 , we have 
L2,N(a 2 ) - L 2 ,N{a\) = (n Ajv) 



1/2 



rJV/ai] 



^ |^(a 2 )r 1 (e 2 (a 2 ,fc)-d 2 (a 2 ,A : )) 

fc=[rAT/a 2 ] 
[(l-r)7V/a 2 ] 

+ ^ [pjv (as) |" 1 (e 2 (02 , fc) - d 2 (a 2 , fc)) - \D N (oi) f 1 (e 2 (ai , fe) - d 2 (ai ,*;))) 
fc=[rJV/oi] 
[(l-r)JV/a 2 ] 

+ Yl \DNiaiT 1 (e 2 {a u k)-d 2 {a u k)) 

k=[(l-r)N/ ai ] 

Then, one remarks that for any finite family / of random variables (Xi) ieI with finite variance we have 

E wl Xl ) = E E i x i x j) < \J EX ? x \J EX j = wl \/ EX ? ) which combined with 

Kiel ) (i,j)ei 2 (i,j)ei 2 Vie/ / 

(x + y + zf < 3 {x 2 + y 2 + z 2 ) implies 



where 



Si = 



E \L 2tN (a 2 )-L 2 . N ( ai )\ 2 < C (NA N )x (Sf + Sl + Sl) 



[rN/ax] 

|^(a 2 )r 1 \ E (e 2 (a 2 ,fc)-d 2 (a 2 ,fc)) 2 

k=[rN/a 2 ] 
[(l-r)N/a 2 ] 



l^vMP 1 J2 \i E (e 2 (a u k)-d^ ai ,k)) 2 , 

k=[(l-r)N/ ai ] 

[(l-r)JV/a 2 ] 



S* = 



E - E 

k=[rN/ ai ] 



I^MP 1 (e 2 (a 2 , k) - (P(a 2 ,k)) - I^Mf' (e 2 ( ai , fc) - d 2 (a u fc)) 



From (jTSJ) , we get e 2 (a, fc) — d 2 (a, fc) = £ 2 (a, fc) + 2e(a, fc) d(a, fc). Moreover, the random variables X — e(a, k) or 
X = d(a, fc) are centred Gaussian random variables, thus we have y/E (X 4 ) = V3E (X 2 ). Then, by combining 
this remark with Cauchy-Schwarz inequality and Lemma lA. 21 we deduce 

E (e 2 (a,k)-d 2 (a,k)) 2 < C {E e 4 (a,fc) + E e 2 (a,k) d 2 (a,k)} 

< C |se 4 (a,fc) + V / Es 4 (a,k) x y/Ed 4 (a,k)} 

< C ^Ee 2 (a,k)y + Ee 2 (a,k) x Ed 2 {a,k) 

< C x <p(N) x |zi(o) + <p(N)Y 

Afterwards 

(N A N ) S 2 < C (NA N ) x <p(N) x |zi(a) + p(N)\ x \D N {a 2 )l 2 x {[riV/aJ - [rA/a 2 ]} 2 

< C (iVAjv) x p(JV) x |a 2 - ai| 2 . (53) 

since \D]y(a 2 ) \ ~ (1 — 2r)N/a as A~ goes to oo and Ii(a) is bounded. The same calculations provide 

{NA N )S% < C {N A N ) x tp(N) x \a 2 - a x f . (54) 

Next, we derive the upper bound for S 2 . Let us stress that the functions a i— > 2V |Z3Ar(a)| _1 converges uniformly to 
a(l — 2r)~ 1 when 2V goes to oo. Thus one can replace Sf by 5| where 



l(l-r)N/a 2 ] 



S* = 



(1 - 2r)A^ 



\JE[a 2 {e>{a 2 ,k)-d?(a 2 ,k))-a L (e 2 (ai,k) — d 2 (ai, fc))f 



fe=[riV/ai] 
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Then, by using (|15p , we get the following expansion of the term define below 



fk ■= a 2 (e 2 (a 2 ,k) - d 2 (a 2 , k)) - a\ (e 2 (a 1 , k) - d 2 {ai 1 k)) 

= (V^2 e(a 2 , k) - y/al e(oi, fc)J x (^y/a^ e(a 2 , k) + ^/a^e(ai, k)\ 

+2 y/a2 d(a 2 , k) x (^^/a^ e{a 2 , k) — v / a7e(ai, fe)J + y'oi e(oi, fc) x ^^"2 d(a 2 , fc) — v^i d(ai, 

We lay the emphasize on the fact that all the random variables in the above formula are Gaussian centred variables. 
But for two Gaussian centred random variables, say X and Y, we get E (X 2 Y~ 2 ) < \JE X 4 x \JEY 4 = 3 {EX 2 ) x 
(EY 2 ). By combining this remark with Lemma TA. 21 one obtains 

Efl < C i^E (ya^e(a 2) k) - v / oIe(a 1 ,fe)^ x E [yE 2 's{a 2 , k) + y/a{ e(a 1; k)j 

+ a 2 E d 2 (a 2 , k) x E (^/a 2 e(a 2 , k) — y/aie{ai, k)j + ai Ee 2 (ai, k) x E (^/a 2 d(a 2 , k) — y'oi d(oi, k)j | 
< Ctp{N) \a 2 - ai| 2 {aiXi(ai) + (o x + a 2 ) p(JV)} + C <p{N) E (y^d(a 2 ,k) - ^ d(oi , fc) V 
But Taylor Formula implies the existence of a real numbers 9t € (ai, 02) such that 

v / a^d(a2,fe) - •v/ai^ffli)*) = (02-01) J (jjP) ^ (Jf ~ kA r^j Xlt) dt 

and after E (^/a^ d(a 2 , k) — y^oi d(<xi, fe)^ < C |o2 — ai | 2 . Indeed, one observe that since t 6 (a m i n , a max ), one 
haves l/# 2 < l/a 2 nin . This implies 



£ ( y y/a 2 'd{a 2 , k) - y /ald{a 1 ,k)\ 



< 



\a 2 - 01 



a a - ai| 



ip' Q - fcA w ) V' (| - fcAw) 5(0, u, 0, «) du dv 



(|u| 1 +**l, t < 1 + |u| 1 +*') du) 







/ sup 


^'Q-kA N ) 


V/iR e€(o mm ,a m „) 





On the other hand, the fast decreasing of the function tjj' insures 

\u\ 1+Hk 1 u <i + \u\ 1+Ho ) du< 00. 







/ sup 




/ JRee(o m i„,a max ) 





Therefore, since ai, a 2 , Ti(a) are bounded and <^(iV) — * as N goes to 00, we have E f 2 < C (p(N) \a 2 — ai\" 
This leads to 



([(1 - r)N/a 2 ] 



rN/ax 



(NA N )-S 2 < C N~ ■ (N A N ) x (p(N) \a 2 — ai[ 
< C ■ (N A N ) x ip(N) \a 2 - ai \ 2 . 

Eventually, combined with (l53l [54]) , one obtains 

E \L 2:N (a 2 ) - L 2 , N ( ai )\ 2 < C ■ (N A N ) x <p(N) |a 2 -ai| 2 . 

But, Lemma [A.2I iii) implies that (N An) x (f(N) converges to when N converges to 00, therefore we deduce (|52|) . 
This finishes the proof of the tightness of the sequence (Ljv(a)) a . < a<a ■ Now, the functional Delta method 
(see for instance Van der Vaart, chapter 20, p. 297), provide a central limit theorem for log(//v(.)) — log(Ii(.)), 
because the function log(.) is a Hadamard-diffcrentiable function on the space of cad-lag function on [a m i n , a max ]; 
this completes the proof of Theorem 13.11 ■ 
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A. 2 Proofs of section 4 



Proof. [Proposition 14. lj We lay the emphasize on the fact that, in this proof, we generalize the 
choice of the frequencies by considering a^r = (NA;y) q , with q > 0. 

For a given N, denote T* = = 0, t*, ■ ■ ■ , t* K , t* K+1 = a N ) such as : 

/«< — </«+!, for all j = l,...,K 

3 a 3 

and for T=(0, ti,--- ,t K ,a N ) £Ap\ we denote z\ N) = s/NA N <y t - log Ji(l// t -)), 

Y ]tj,t j+ i] = i Y t s +i,--- , Y t s+ i-T N y, x ]tj,t j+1 ] = Qoefy+%,i)i<i<(t :i+1 -t j ), z ]t^l J+1 ] = ( z tf+ir-- , z tf + \-T N Y- 

First step : We would like to prove : cDj; -—> lo* for all j = 1, ... , K. 

Denote = Q^(t*,A(t*)) where A(t*) is obtained from a linear regression of (Y^) on (log/j) for i — 

t* + !,-■■ ,t*+i - tn. Let e > and II T - T 1 ||oo= max It,- - t'A for T = (0,ti, • • • ,t K ,a N ) £ A^p and 

je{i,-,iC} 

T' = (0,ti,.-- ,t' K ,a N ) & A ( p. Then, we get, 

p(\\T-t* ||oo>eajv) <P^min QW^A^^Q^ 

where V eajv = |t € -4^°, || T - 1* \\oo> eOjv}- We want to show that for all T £ F eaj 
= o(QW(T,A(T))). In fact, 



A'+l 



NA 



N 



3=0 
A'+l 



7 W 
Ht*,t* ,i 



< 1 y f7 W y 7 W 



1 



[l,a N ] ' 



From Proposition 13. 1[ we deduce 



13 (af. 



aN 1 J L J JV— oo Jo \Jmin \PJn 



du. 



(55) 



which is a positive and random variable because Z is a continuous Gaussian process. Afterward, for a sequence 
{"4>k)k S and a sequence of random variables (£fc)fceiv, we will write £jy = Op(V'iv) as TV — > oo, if for all £ > 0, 
there exists c > 0, such as , 

p(\€n\ <c-^n) >l-e, 
for all sufficiently large AT. Here, we obtain : 



\ NA 



N 



(56) 



Now, let T £ V eaN , we want a lower bound of QW(T, A(T)). Wc use the following decomposition 

XjAj — log 2^(1/ /., 
2 [Fi - logli (!//„)] x [A, A, - log2i(l//() 



K + l tj + l—TN 

QW(T,A(T)) = E E K-fogii(i//i)] : 

j=0 i=tj + l 



^3- 
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Then : 



K+l 



N 



1. Since Q t = ^ as P reviousl y we § et 

'/3/m 



NA 



N 



(57) 



2. Let r = log 



min < log 

3=1,-, K\ 



Then, for all j € {0, 1, • • • , if}, - > i* + ra N . 

Since T G V eaN , we have 77 = min{e, r, log(/3/a)} > and there exists an integer j £ {0, ■ • ■ , if + 1} for which 
there are no estimated abrupt change in the interval [t* — rjciN, t*] or [t* — tn ,t* — tn + na^]- Thus there 
exists fee {0, • • • , K + 1} satisfying [£* — r/aN, t*] C [</c, tk+i — tn] (we follow here a similar proof than Bai 
and Perron in Lemma 2, p 69) and 

Q 2 > I^Afc-logX^l//,)! 2 

i=t*-ija N + l 



i=t*— »?a N + l 



I 

a N 



(58) 



with : 

• A(ff,cr) = log (a 2 • tf ff (V>) ) - (2// + J 1 • 1..- 

• B(H,a) = -(2#+l)-log 



/3 



for all (if, cr) G K; 



Pfrr 



log 2i(l//<) -log Xi 



for all (H,a) £ JC; 

. frain V ^/mi 



Since for all (ii, cr) G /C, the function x 1— > L(H,a) i x ) = (.MH, a) + x ■ B(H, a) — g(x)\ is an infinitely dif- 
ferentiable function on M, we know from the theory of Riemann sums that : 



u N {H,a) = — V 

& N ' 

i=t* -r;a N + l 



A(H,<j) + — -B(H,a)-g( — 



N- 



u(H,a) 



(A(H, a)+x- B(H, a) - g (x)) 2 dx, 



S--7) 



with s* = log 



UJ"i 



f„ 



log 



( Oijr, 



\Pfn 



lim Moreover, the sequence (un(H,<j))n converges uni- 



formly to u{H,a) because for N large enough 

sup \u N (H,a) - u{H,a)\ < ( + V 



a N 



■ sup 



sup 



(H,a)<EK. |0<i<(s^+1) 



dL 



dx 



(x) 



N- 



0. 



dL 



dx 



(x) 



< 00. As a consequence, 



since K, is a compact set of [0, 1] x]0, oof and thus sup < sup 

{H,<r)eK {0<x<(s' K +l) 

from (I58|) and since we assumed that (H^ € K, for alH = 0, • • • , K, for some sufficiently small, fixed 

£ > and for all sufficiently large N, 



Q2 > aN 



A(Sl N) A N) ) + - ■ B(Hl N \*l N) ) -9(x)) dx - £ 



(59) 
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But it is impossible that there exists (a, b) € 1R 2 such as g(x) = a + b ■ x for all x £ [s| — 77, s^], i.e., 
Xi (ci • e C2X ) = e a ■ e b - x for all x £ [s* - 77, s*} with ci = , c 2 = log ( J , which can also be written 

^ ' J vain \ P J max J 

as : 

TiO) = ai • x bl for all x e [a/u^a/w* + rj'}, (60) 
with 77' > and {a\, b\) £ M 2 . Indeed, assume now (j60l is true. But, for all x £ {a / uj* , a / uj* + rf], 

TM of ** ™* i+l l^")l 2 * o_ ^ 2^ + i ^ ^ 

= 2 I aj_ x • x j- 1 ^ / -^jp — — du + cTj • x ^ j 2 h*+i du I . 



Then -—Lta/u*) = a x ■ —^—(a/uj*) for n = 0, 1, which implies that &i = (2ff* + 1) and ai = 2a* 2 K H .Nj) 
ox 11 J dx n J j j 3 

(here, we use the equality ip{a) = 0). Thus, for all x £ [a/uj*, a/uj* + r/], 

a* 2 r 1 1^)1" Ju - r^jSlO!!.^ 

(a* 2 cr* 2 \ 

T7+T - ~^nh+i )dy = 0, 
y y 3 / 

f (7* 2 = (T* 2 

and hence < if* _ ^ . But this condition is impossible from Assumption (Bk) and consequently there 
is no (a, b) £ Ft 2 such as g(x) — a + b ■ x for all x £ [s* — 77, s*}. 

The function g belongs to the Hilbert space JL 2 ([s* — 77, s*];dx). Since C = {A + B ■ x, x £ [s* — 
77, Sj], (A,B) £ M 2 } is a closed linear subspace of K, 2 {[s* — 77, s1j];dx), there exits a distance between g 
and £ in JL 2 ([s* — 77, s*]; dx), i.e. there exists (A, B) £ Ft 2 such as 

[ 3 (A + B-x-g(x)Ydx= inf (A + B ■ x - g(x)) 2 dx = C > 0, 

because g yt C Then, by choosing £ such as < £, < C/2, the inequality (|59[) implies : 

Q 2 > y ' ffi /v ( 61 ) 
for all sufficiently large N, with C a real positive number only depending on 77, s*, H*_ 1 , H*, crjLij a* and 

3. The previous evaluations of Q\ and Q2 provide an upper bound of Q3.We get 

/K+lt k + 1 -r N \ X/ 2 

Q 3 < 2(gi) 1/2 ^ ^ (^Afe-logX^l//,)) 2 

\fc=0 i=t k +l / 

, . 1 I 2 

< 2(Q 1 ) 1/2 x [a N - sup <|2sup{(log/,l)-A) 2 + 21og 2 X 1 (l//) 

\ f min <f<f m ax I XEK 

= Opf-^V (62) 

We deduce from (|57)) . ([6T|) and (|62|) that Qi = o(Qa) and Q3 = 0(^2), which implies 
P ( min g (Ar) (T, ACT)) > — ■ a N ) — > 1 and thus 

lim P f II T — T* IIoq^ sax) = =^ S, (JV) w*. 



TV— *oo V / AT— >oo 
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Second step : For j = 1, ■ ■ ■ , K, we want to prove that if 3/4 < p < 1 and < q < 1, for all e > 0, there exists 



< C < oo such as for sufficiently large N, P (a} N 



~(JV) 



> C < e. 



Mutatis mutandis, we follow the same method as in the proof of the convergence in probability. Now, let < p < 1, 

< 77 = -min{T, log(/3/a)} and consider min TeW v Q( N ~>(T,A(T)) with 
2 c «» 

W^ = {Te^ ) ,Ca^<||T-r ||oo<?7ajv}. 

Then, as previously, for T e W£ aP and AT large enough, it exists j e {1, • • • , K} such as 

tj + Ca p N < t* < t j+ i - t n 

(the following proof is valid even if one considers the alternative t* < tj — Ca p N ). Then 

tj + i — tn 



QW(T,A(T)) > J2 (Y i -logl 1 (l/f i )) 2 + (X i X j -logl 1 (l/f i )) 2 + 

+2(Y t - logJiCV/OJ^Aj - IogJi(l//i)) 



i=t;+i 



1. First, we have again, 



a N 



NA 



N 



tj + l—TN 



(63) 



i=t*+l 



Moreover, for a* = i e {tj + 1, ■ ■ ■ , i?} and N large enough, a* ~ a/u)*, and 



Ji(ai) = Ji — + a< - — J( — + O Ui - — 



(64) 



2. Secondly, Q 2 = (XjAj - log Ji(l//i)) 2 . But wc know log Ji(l//<) - XA* for i e {t* + 1, ■ ■ ■ , t j+ i - t n }. 



But Ji( 0i ) = 2 ^ioJ H; - 1+1 

IogJi(l//i) = X;A* + 



and 



2H? 



; thus for i G {t, + 1, • • • , t*}, 



fmin , / Pfmax 



(2ff* + l)^l g 



/3 °W> 



Oat 



(Zjv 



Then, with Aj = (aj,bj)', one gets for i e {t* + 1, • • • , tj + \ — tat}, 

(X-Xj - logli (!//<)) = (log/i - 1^7) (oj - a*) + Z, 



XXX indicates the empirical mean of XXX between tj + 1 and tj + i — t^. Thus, 



tj + l —T N 

Q' 2 > E ((log/ i -loi/)(a J -a*) + 



»=t;+i 



1 - 



VN^n 



(65) 



(66) 



(67) 
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and thus. 



tj + l— TJV 

(log/.-IoiT)^-^ 7 ) 

tj + l—T N 

£ (log/.-Ioi/) 2 
J2 (log /< - loi/) Iog2i(l//<) + ^=_Z s (Ar) - ToixT 



JV 



Vna n 



-Z 



J] (log/i-loi/) 2 

i=tj+l 



(log /< - log 7) (logTxCV/i) - x;a*) 
^ (log/.-bgT) 2 



ij + i — t~n 

J2 (log/.-bgD^w-z 

1 i=t,-+l 



\/iVA 



iV 



tj+i — tn 

(log/i -bgT)' 

i=t,+l 



From the definition of (log/i), 



£ (log/.-Ioi/) 1 



log 



12 D \4 



/3/n 



(tj+i - t n - tj) = O(ajv). 



Expansions and (p5|) imply there exist two constants Ci > and C*2 > such as for AT large enough 

t* 



t*. — t ■ 
a N 



< 



(log fi - bgT) (logXi - XiX* 



Y (log/i -loi/)' 

i=t,+l 



< C 2 



a N 



Moreover 



y/NA 



tj+1 — TN 

£ (log/.-bgD^W-z 

1 t=t,+l ( \ 
— Up 



AT 



tj+1— TN 



Y (log/, -log/)' 



Thus, we deduce from (l68l) that 
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As a consequence, for (p, q) such as 4g(l— p) < 1 (for instance, p = 3/4 and q = 1), then ( — | • \/ NAn > C , 

V ffljv J 

and thus for all e > 0, for iV sufficiently large, we can chose C > such as : 

P ^ (log/,; - fog7) 2 (^p) < ((log/, - fog?) (%• - a*))^ > 1 - E. (70) 

Now, from (O, (TTOl) and with P(< m - r w - <* > OA r) — > f , for (p, g) e [3/4, 1] x [0, 1], for all e > 0, for TV 
sufficiently large, we can also chose C > such as : 

\ aN / %=t*+i j 

=> P (C 4 • C 2 ■ a%- 3 < Q' 2 ) > 1-e, (71) 
with C2 > a real number not depending on C, N and e. 

3. Finally, from the classical bound of Q' 3 , we obtain, 

Q , 3 <2-(Q' 2 ) 1/2 -(Q' 1 ) 1/2 . 

But, following a similar method as previously, from ([701 one can find a upper-bound for Q' 2 , i.e. for (p,q) £ 
[3/4, 1] x [0, 1], for all e > 0, for TV sufficiently large, we can also chose C > such as : 



with C3 > a real number not depending on C, TV and e. Thus, for (p, q) £ [3/4, 1] x [0, 1], for all s > 0, we can 
also chose C > such as : 

(2p-2 \ 
^'^■tM - '-■ (72) 

with C4 > a real number not depending on C and N. 

Now, from ([64]) , ([71]) and ((72J), one deduces that for (p,q) £ [3/4,1] x [0,1], for all e > 0, for TV sufficiently 
large, we can chose C > sufficiently large such as : 

P [ min QW(T, A(T)) > C 4 • ^ • a^ 3 ) > 1 - e. 
and thus like Ql W) = P from (J5BJ), 



that leads to P 



( a N P 



~(JV) * 



P I min QW(T,X(T)) < Qi N) | < e, 



> C J < e for sufficiently large C and TV. 



Proof. [Proposition 1472] From Proposition 14. 11 we deduce that Vj = 0, • • • , TV, 



([f,f]c[%r^]) 



— > 1. 

JV^oo 
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Denote A\ N) the event [tjf\V^ N) ] C [t* p t* j+1 - t n }. Then, Vj = 0, • • • , K and V(ss, y) € -K 2 
F (V^VA^ (Xf ' - A*) e] - oo, s] x] - oo, yf 



(Af ! ) x P (v^VA^(Af > - A*) e] -oo,x]x]- oo,y] \ A^) + 



+ P x P (^NA^ (Af - A*) e] - oo, s] X] - oo, y] | Ap) . 

Now, since P (^/WK^ (lf ] - A*) G] - oo, x] x] - oo, y] \ A^ N) ) < 1 and P (A^ N) ) =1-P (^f ) > we obtain : 
P (Af >) • P ( v^VA^ (Af } - A*) e] - oo, x] x] - oo, y] | Af >) 



< P ( V^Ajv ( Af > - A*) g] - oo, a] x] - oo,y] 



< P 

s(jV) p 

— r ^ 

OO 

^ ( fl *(fc))i<*< m aH 
orem JT7J), for all (^fe)i<fe< m € iR m , we get 



(v^VA^(aW-A*) g]- oo, a:]x] -oo,!,] I if'j+l-Pff). (73) 



Since (2>- — ► and — ► ; therefore (tD,- — > (uj* ,u}*,-,), we have 

J N^oo J J N^oo J J J N->oo J J 

ifk)u^iTTW \A N )\ (Sit^lKKm and X^' X*. Thus, from Proposition |3~T1 and central limit the- 



/ m \ / m 

P v^A^ (f} N) - Xfh*) ef[] oo, x k ] | AW P U- G f[] - oo, x k ] | A\ 



k=l / \ k=l 



— ► o, 



with 2 A/" m (0, S*) and E* = ( cov ( Z (^77t) > Z (-^7T )) ) ( it; explains the expression (J24]) of E*). 

From the equality Af° = ((Jfj^J'xj" 5 ) (x!j N) )'Y} N) , we deduce that for all (a;, y) € JR 2 , with ~ JV 2 {0,T^) 

and if = (a; a;) \y;v;a; (a; a;) '. 

P (a?° - A*) G]-oo, x]x]- ex., y] \ A { f ] )-P e] - oc,z]x]-oc, y] | if) — * 0. (74) 



We also have 



(li e]-oo,a:]x]-oo,y]) + P (if } ) 



- 1 < 



/_ (m \ p {lj G]-oo,a;]x]-oo,y]) 

<P(^G]-oo,x]x]-oo,y] \Af ) ) <^ — — ^ (75) 

V ' P [ A j ) 



Now, as P (a\ n) ) — > 1, from (FFl, J71 and (f75j), we deduce that for all (x, y) E R 2 : 

V J J N^oo 

P (y/NA N (Xf > - A*) G] - oo, s]x] - oo, y]) — + P (|, G] - oo, x] x] - oo, y]) , 



that achieves the proof. 



Proof. [Proposition [473] First, from the expression of each s k i given in ([24]) and with Ai m (M) the set of real 
m-by-m matrix, the function E : (H,u,v) i— > T,(H,u,v) G M m {IR) is a continuous (and therefore measurable) 
function of (iJ, u, u) for in a compact set included in ]0, 1[ and (u, w) &]f m in> fmax[ 2 - For all j = 0, • • • , X, we 
have : 



34 



Statistic of multi-scale fractional Brownian motion 



1. from Assumptions (B K ) and (C), (h/ N) , erf ) € /C and (D J (JV) , ^+1) e]/mm, /maip! 

2. from dH]) and flU, h/ N) Hf, wf lu*, wf and therefore 

N-^oci iV— *oo TV— >oo 

(JV) ^(jv) ~(AT)n P /XT* * * N 
N— >oo 

As a consequence, — ^{Hj ,Oj N \Qj+i) S|, for all j = 0, • • • , K, and since S(.ff, u, v) is an invertible 

TV— >oo 

covariance matrix for all (H 7 u, v) £]0, l[x]/ m ^ n , / maa :[ 2 , 

(£f } ) ~* (s*) _1 , for all j = 0, - ■■ , K (76) 



Secondly, denote 



Mf ) = ((If )y( S |) xf ^"'(xf >)'(e£" ' 



The 2-by-m matrix M"f ^ verifies : 

Af = W = a; + _2=Mf >zf > 

with if = (zWh/fA) and Z,- = fz(l/p|(fe)) N ) from the central limit theo- 

rem ([XT]) . In the same way, 

Af > = M^Y} N) = A* + -J=jgf 



3 ^fWK? 3 3 



From (7BJ), we obtain m\ N) - m\ N) 0, and thus, 

3 3 N^oo 



V J J I J J iV-»oc 



with Aff ^zf ^ — 7^(0, T 2 J ) (the same covariance matrix as that obtained with a generalized least squares 



estimation), and this implies Proposition 14.3 



Proof. [Proposition 14.41 For each j = 0, • • ■ , K , one first show that 



na n .\\yW-xW A f)||| w ^ xV _ 2) . (77) 



Indeed, || Yf > - If A f ||| w =|| ll| ( ^ ^ II ^f > ll| (N) where pf > = 7 m - xf >^> is 

the matrix of the orthogonal projector in lR m on the orthogonal of Vj, where Vj = {xf A, A S -S? 2 } is the 
2-dimensional subspace of lR m generated by xf (here the notion of orthogonality is based on the inner product 

< it, v >g(N)= u' ■ (sf ^ ■ t) for ii, u £ _ZR m ). From the previous proofs, we know : 
. £f) s* if) _^ X* and therefore pf P* ± where 

3 N^oo 3 3 N^cc 3 3 jV->oo 3 

-1 



P? X = I I, 



X* (x* (£*) 1 X J *^ X* (£*) 1 J is the matrix of an orthogonal projector on a (ra - 2)- 
dimensional subspace of M m ; 



v 

<u,v >s(«) — > < u, v >s* for u,v G J? m ; 
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• z\ N) with Zj 2 A/" m (0,E*). 

Consequently, || P^Z-* ^> || P/ X Z,- |||*. From Cochran's Theorem, we know || P/ ± Zj |||, 

~ X 2 (m — 2) and therefore (|77|) is proved. 

Moreover, with the notations of Proposition 13. 1[ if log/ > log/' + log/3/a then cov(Z(l//), Z(l//')) = 0. But 
for all S {0, • • • , K} 2 , i ^ j, Vfc G {£/f°,- • • , V} N) } and Vfc' e {t/f , • • • , vf°}, | log/ fc - Iog/ fc /| > log/3/a. 
Thus, we deduce that the different Aj- are asymptotically Gaussian and independent. It provides the end of the 
proof of the Proposition 14.41 ■ 
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